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An estimate of the net direction of climate interactions in different geographical regions is made 
by constructing a directed climate network from a regular latitude-longitude grid of nodes, using 
a directionality index (DI) based on conditional mutual information. Two datasets of surface air 
temperature anomalies—one monthly-averaged and another daily-aver aged—are analyzed and com¬ 
pared. The network links are interpreted in terms of known atmospheric tropical and extra-tropical 
variability patterns. Specific and relevant geographical regions are selected, the net direction of 
propagation of the atmospheric patterns is analyzed and the direction of the inferred links is vali¬ 
dated by recovering some well-known climate variability structures. These patterns are found to be 
acting at various time-scales, such as atmospheric waves in the extratropics or longer range events 
in the tropics. This analysis demonstrates the capability of the DI measure to infer the net direc¬ 
tion of climate interactions and may contribute to improve the present understanding of climate 
phenomena and climate predictability. The work presented here also stands out as an application 
of advanced tools to the analysis of empirical, real-world data. 
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Information-theoretic tools are used 
to construct directed climate networks 
from time-series analysis of observed 
climatological data. Specifically, sur¬ 
face air temperature (SAT) anomalies 
are considered. Two datasets—one 
monthly-averaged and another daily- 
averaged— are used. Directed links 
among the network nodes are defined 
via an analysis of the net direction of 
information transfer. A predictability 
measure—based on conditional mutual 
information—quantifying the amount 
of information in a time-series x(£), con¬ 
tained in r time units in the past of an¬ 
other time series, y(t ) is used. The re¬ 
sulting directed network is then studied 
and a full agreement with state-of-the- 
art knowledge in climate phenomena 
has been found, validating this method¬ 
ology for inferring the net directionality 
of climate interactions, directly from 
the data. No weather assumptions or 
models are made, except for the appro¬ 
priate setting of the parameter r which 
is sensible to the shorter or longer auto¬ 
correlation of the time series. 


* juan.ignacio.deza@upc.edu. 


I. INTRODUCTION 


Network theory is a well-known framework for describ¬ 
ing complex systems composed of many interacting com¬ 
ponents EH1- Many systems can be straightforwardly 
represented in terms of a well-defined set of nodes cou¬ 
pled among them via links that have a clear physical 
interpretation. This is the case, for example, of airport 
networks [E|, of social interactions |6j or the Internet [7|, 
just to name some. In other systems it is not clear how 
to define the relevant nodes, and/or there is no obvious 
interaction that can be used to define links. An example 
of this situation is the Earth climate system, in which 
the lattice of grid points from measurements or models, 
is defined to be the network set of nodes, and it depends 
on the resolution of the dataset analyzed. Many climato¬ 
logical fields—such as surface air temperature (SAT) or 
the geopotential height (GH) at a certain pressure level— 
can be used to define links via an analysis of significant 
correlations |M3- 


Climate networks (CNs), which represent the statis¬ 
tical similarity structure of spatio-temporal resolved cli¬ 
matological variables, depend on the definition of nodes 
and links 0. Their regular spatial sampling results in 
a small-world topology [12j, and thus, a careful interpre¬ 
tation of the inferred network is required. Nevertheless, 
CNs have been successfully employed to analyze climate 
features including the global connectivity jl3l4l5l |. the 
identification of community structures [16|, [l7| and the 
study of the possible collapse of the meridional overturn¬ 
ing circulation (MOC) 
others. 


18| on the north Atlantic, among 


A relevant drawback of this correlation analysis, 
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which uses symmetric similarity measures (such as cross¬ 
correlation or mutual information), is that it yields non- 
directed networks where the presence of a link reveals 
inter-dependency but the direction associated (if any) of 
the underlying interaction is not established. For im¬ 
proving the understanding of climate phenomena and its 
predictability, it is of foremost importance not only to be 
able to infer the presence of a link between two nodes, 
but also, to infer the direction of this interaction. 

A path to overcome this limitation is by constructing 
weighted climate networks, where the weight of each link 
is composed of two numbers: the correlation strength 
between the two nodes and the time delay which max¬ 
imizes this correlation strength. The sign of the time 
delay gives information about the direction of the link. 
Using this approach, the high sensitivity of the network 
links to El-Nino events was demonstrated, even in geo¬ 
graphical regions far from the Pacific ocean [3 Hi 

An alternative approach for assessing the directional¬ 
ity of climate interactions involves the use of Granger 
Causality. In (2l| it was argued that the inclusion of 
several variability patterns like NAO (North Atlantic 
Oscillation), PDO (Pacific Decadal Oscillation), ENSO 
(El Nino-Southern Oscillation), and NPI (North Pacific 
Index)—which occur naturally and explain an important 
part of the global atmospheric variability—into a nonlin¬ 
ear network-like prediction method, largely improves the 
predictability of global temperature over seasonal time 
scales. This is so even when the time scales of the pat¬ 
terns used of in the order of decades. The study sug¬ 
gested a causal directional influence among these major 
oceanic and atmospheric modes and global temperature 
variability, not only over their own time-scales but also 
over much shorter ones. Granger causality has also been 
used to test interdependence among ENSO and the In¬ 
dian monsoon [22|. A non-symmetric bidirectional and 
even alternating character of coupling was found that ex¬ 
tends previous knowledge about the presence of negative 
correlation and intervals of phase synchronicity between 
the processes. 

A third approach for directionality detection is based 
on information-theoretic measures [23|. For example, 
in [24| information transfer from larger to smaller time 
scales was detected in daily-averaged SAT time series 
as the causal influence of the phase of slow oscilla¬ 
tory phenomena with periods of about 6-11 years in the 
amplitudes of the variability, characterized by smaller 
timescales, from a few months to 4-5 years. 

The reliability of directed climate networks detected by 
bivariate nonlinear methods based on information theory, 
compared to those generated by linear Granger causality 
analysis, was studied in [25j. Several algorithms for esti¬ 
mating transfer entropy with a wide range of parameter 
choices were considered. As transfer entropy is a special 
case of the conditional mutual information, it reduces to 
Granger causality for linear Gaussian processes, and it 
usually requires longer time-series for accurate estima¬ 
tion. It was shown that all the causality methods consid¬ 


ered provided reproducible estimates of climate directed 
interactions, with the linear method outperforming the 
nonlinear ones in terms of reliability. 

In the present paper, directed climate networks are 
constructed by using a predictability measure which, to 
the best of our knowledge, has not yet been applied to cli¬ 
matological data. This bivariate analysis quantifies the 
amount of information contained in a time-series, x(t), 
about the past of another time series, y(t ), r time units 
before. More specifically, the directionality index (DI) 
is used, which quantifies the net direction of informa¬ 
tion flow between two time series, and has been success¬ 
fully employed for the analysis of numerical data gen¬ 
erated from coupled oscillators, and empirical data from 
cardiorespiratory recordings and electroencephalographic 
recordings, to name just a few [26l-[28j. 

The objective of this work is to demonstrate that the 
DI methodology for inferring the net direction of in¬ 
formation transfer indeed unveils climatologically rele¬ 
vant phenomena. Clear patterns of variability are un¬ 
covered, both in the tropics and in the extra-tropical 
regions, and their direction of propagation is shown to 
be in good agreement with the current understanding 
of main climate phenomena. The method is shown to 
work with both monthly- and daily-averaged data, and 
furthermore, well-defined atmospheric patterns are un¬ 
covered on daily-averaged data, which are not seen with 
monthly-averaged data. 

The paper is organized as follows. Section lU presents 
the data analyzed and the method used for constructing 
directed climate networks. The statistical significance 
test is also discussed. Section inn presents the results 
and provide a comparison between directed and non- 
directed climate networks. Relevant patterns of global 
atmospheric variability are interpreted from the maps ob¬ 
tained through the DI analysis. Finally Sec. ITVl presents 
the conclusions. 


II. DATA AND METHODS 
A. Data 

We consider two datasets, both corresponding to SAT 
anomalies from the reanalysis of the National Center 
for Environmental Prediction/National Center for Atmo¬ 
spheric Research, (NCEP/NCAR) p]. The data covers a 
regular grid over the Earth’s surface with latitudinal and 
longitudinal resolution of 2.5°, resulting in N = 10226 
grid points or network nodes. 

The first dataset corresponds to monthly-averaged 
SAT data. Since the data cover the period from January 
1949 to December 2013, in each node there is a time series 
of 780 data points. The SAT anomalies are calculated as 
the actual temperature values minus the monthly aver¬ 
age, and they are normalized by the standard deviation. 
Each time-series is linearly detrended. 

The second dataset consists of daily-averaged SAT 
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from the same source, with the same spatial resolution, 
and covering the same period of time (the length of each 
time series is 23725 data points). In this case, for cal¬ 
culating the anomalies, the daily average has been sub¬ 
tracted from the actual temperature value, and the leap 
days have been discarded. The time series are also de¬ 
trended and normalized by the standard deviation. 


B. Directionality measure 

We consider the SAT anomalies time series in two 
nodes, X(t) and T(t), which are characterized by proba¬ 
bility distribution functions (PDFs) px, Py, and by their 
joint PDF, pxy- To assess the directionality of the link 
between these two nodes the directionality index (DI) as 
defined in [2jJ 30| is used: 


DI X y(r) 


Ixy{t) - Iyxjr) 
Ixy ij ) + Iyx{t) ’ 


(1) 


where Ixy(t), Iyx{t) (conditional mutual information) 
are defined as: 


Ixy(t)=I(X;Y \X t ) 

= H(X\X t ) + H(Y\X t )-H(X,Y\X t )- (2) 
Iyx(t) = I(Y;X\Y t ) 

= H(Y\Y t ) + H(X\Y r ) - H(Y,X\Y r ), (3) 

with X T = X(t - r), Y r = YU- r) and H(X\Y) being 
the conditional entropy |26|, |30| . 

In terms of three time series X(£), Y(t) and Z(£), 
I(X;Y\Z) indicates the amount of information shared 
between X(t) and Y(t), given the effect of Z(t) over 

Y ( t ). In this way, to assess the information transfer from 
X(t) to Y(t), Z(t) is replaced by the past of time-series 
X(t). Thus, I(X;Y\X r ) quantifies the amount of infor¬ 
mation shared between X(t) and T(t), given the influ¬ 
ence of X(t — r) over Y(t). Analogously, to assess the 
information transfer from Y to X, Z(t) is replaced by 
the past of Y(t). This special case of conditional mutual 
information is also referred as transfer entropy. 

The directionality index, Dixy, then quantifies the net 
information flow. From the definition of Dixy, Eq* (HJ) , 
it is clear that Dixy = — Dlyx- Also, — 1 < Dixy < 1: 
Dixy = 1 if and only if Ixy 7^ 0, Iyx = 0 (z.e., the 
information flow is X —>• Y and there is no back coupling 

Y X) and Dixy = — 1 if and only if Ixy = 0, Iyx ^ 0 
(z.e., the information flow is Y X and there is no back 
coupling X —>Y). 

Naturally, r > 0 is a parameter that has to be tuned 
appropriately to the time-scales relevant to the particular 
dataset. If r is too small DIx y(t) will capture short 
time scale directionality, and may fail if the time series 
behave too similarly on those time scales as they do if 
they are subjected to the same external forcing. On the 
other hand, if r is too large, larger than the decorrelation 
time of the time-series, the effect of the past X over Y 


(and of Y over X) will be negligible and DIxy(t) will 
be a small and in principle random value. In the next 
section the criterion employed for accessing the statistical 
significance of the DI values will be discussed, and in 
section ED a thorough study of the effect of varying the 
parameter r will be performed. It is worth noticing that 
in the two datasets considered, r = 1 has a different 
meaning: in the monthly-averaged SAT time-series, the 
minimum value of r is one month, while in the daily- 
averaged data, the minimum value of r is one day. 

To compute the PDFs associated to each time series, 
and the joint PDFs a 10-bin histograms of values is used. 
An alternative approach is to use a symbolic transforma¬ 
tion known as ordinal analysis [26|, l3ll l33l | . A preliminary 
study suggests that using ordinal analysis offers as an 
additional advantage the possibility of finding the direc¬ 
tionality of the links at different time-scales; this study 
is still in progress. 

It is important to note a drawback of using the direc¬ 
tionality index for network construction: it does not dis¬ 
tinguish indirect from direct information transfer. Given 
three time-series, X, Y and Z, positive and significant DI 
values Ixy > 0 and Ixz > 0 might also give a significant 
value of Iyz —either positive or negative—even if there 
is no “direct” information in Y about the future or the 
past of Z, as the information will have been “indirectly” 
contained in X. 


III. RESULTS 

Firstly, the results of the statistical significance anal¬ 
ysis of link directionality, when using monthly-averaged 
SAT datasets, are presented. Monthly datasets are used 
in order to to compare with the undirected networks re¬ 
ported in [33l. Afterwards, the results of the analysis 
of the daily-averaged datasets are presented. They are 
consistent with the monthly datasets results, furthermore 
uncovering additional patterns of atmospheric variability, 
not observed in monthly data as the averaging procedure 
filters high frequency and variability. 


A. Statistical significance analysis 

To address the significance of DI values, 100 surro¬ 
gates were generated using the bootstrap (BS) algorithm 
[34{ . The null hypothesis considered for MI is that the 
processes are independent from each other. The null hy¬ 
pothesis for DI is that there is no preferential direction 
of the interaction. 

The BS algorithm randomly resamples with replace¬ 
ment from the original datasets using blocks of data of 
approximately the size of the autocorrelation time of the 
time series, and then computes the estimators (MI and 
DI) from the resampled data. Doing so, both the statis¬ 
tics (histogram) and the power spectrum of the original 
time series are approximately preserved. In this way, 




4 




FIG. 1. (Color Online). Procedure of constructing significant directionality maps from raw DI calculations. In panels (a) and 
(d) the unfiltered DI maps are shown for two nodes, one in the central Pacific ocean and one in the Indian ocean, indicated 
with triangles. In panels (b) and (e) only the statistically significant MI values are shown. These results are combined in (c) 
and (f) where only the links that have both, statistically significant MI and DI values are shown. 


two different empirical distributions (one-tailed for MI 
and two-tailed for DI), were obtained for each link, and 
significance thresholds were extracted from them. 

Afterwards, for each link, the MI value was calculated 
from the original datasets and a significance test was ap¬ 
plied. A MI value is considered significant if: MI > ©m/? 
where ©m/ is the threshold derived from the bootstrap 
MI distribution. In addition—for the MI significant links 
only—the DI value computed from the original dataset 
was compared to the DI thresholds ©J> 7 and @^j. The 
DI value was considered significant if DI > or if 

DI < @^j. This two-step significance test was performed 
with the scope of assessing firstly the presence of a link, 
and afterwards its directionality. 

A graphical explanation of the full procedure is shown 
in Fig. [l] where, for two nodes (one over the Pacific and 
one over the Indian ocean) the unfiltered DI maps are 
displayed in panels (a,d), the significant MI values are 
displayed in panels (b,e) and finally the DI significant 
values for the MI significant links are displayed in panels 
(c,f). The DI maps show positive DI values in red, which 
mean outgoing links, while the incoming links are shown 
in blue. 

The significance thresholds, ©m/ 5 ®di and ©^j, 
extracted from the BS surrogates, were computed as 
@m/ = Umi + 3 x (Jmi, ©£>/ = I^di + 3 x (JDi , and 
= ii d i ~ 3 x <jdi, with fi and a being the mean 
and standard deviation of the corresponding MI and DI 
bootstrap distributions. The process is illustrated in Fig. 


[21 panels (a) and (b). It can be noticed that these dis¬ 
tributions are bell-shaped, which justifies the use of clas¬ 
sical statistics. Also notice from panel [2fb) that for a 
sufficiently large number of BS surrogates, the distribu¬ 
tion tends to be centered around a rather small DI value, 
which motivates the null hypothesis of no preferential di¬ 
rection of the interaction. 

The more time and memory consuming 99nth per¬ 
centile test was also used, and yielded similar results. 

In order to visualize the results of the significance cri¬ 
terion employed, in Fig. [2] panels (c) and (d) display—as 
an example—the DI vs MI values of all the links of the 
two nodes, one located in the central Pacific and the other 
in the Indian ocean (the nodes are the same as in Fig. 
[Tj) . The black dots indicate the disregarded links (either 
because MI is not significant, or because MI is signifi¬ 
cant but DI is not). Significant links are plotted in red 
(outgoing links, DI > ©^ 7 ) and in blue (incoming links, 

di < e~ DI ). 

In Fig. [2j panels (c) and (d) also show that high MI 
values do not imply high DI values. In Fig. [2] (d) one can 
notice that most of the blue dots are located in a narrow 
range of MI values while they are distributed in terms 
of DI values. An inspection of panel (f) in Fig. [1] shows 
that the blue links come to the node in the Indian Ocean 
from a well-defined region in the central Pacific Ocean. 
On the other hand, one can notice in Fig. [2] (d) that the 
few red dots are more distributed in the MI, DI plane, 
and in Fig. [T| (f) the red outgoing links connect the node 
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FIG. 2. (Color Online), (a), (b) Histograms of the MI and 
DI values (t = 1 month) computed from 100 BS surrogates, 
the dashed vertical lines indicate the significance thresholds 
and the solid lines (red online) indicate the MI and DI values 
computed from the original data, (c), (d) Plot of DI (r m 
1 month) vs. MI for all the links of the two nodes in the 
Pacific and Indian oceans considered in Fig. \l\ Incoming links 
are indicated in blue and outgoing links, in red. Significant 
links are plotted in red or blue dots; the black dots indicate 
the disregarded links (either because MI is not statistically 
significant, or because MI is significant but DI is not). The 
triangles indicate the (MI, DI) values of the particular link 
between the two nodes. 


in the Indian Ocean to various regions on Earth. 


B. Analysis of monthly-averaged SAT anomalies 

1. Influence of the parameter r 


FIG. 3. (Color Online). Effect of r on tropical areas using 
monthly-averaged data. In this case a point in central pacific 
(the same point used in the top row of Fig. [T| and in [33|] is 
considered. The values of r are: (a) 1 month, (b) 3 months, 
(c) 6 months, (d) 12 months. Notice the decorrelation of 
the time series for large r. Incoming links are in blue while 
outgoing links are in red. 



As stated in the methodology the correct choice of the 
value of r is necessary for obtaining consistent results. As 
r can be only an integer, using monthly averaged data, 
its minimum value will be of one month. In the trop¬ 
ical areas the influence of the ocean on the surface air 
temperature is a dominant characteristic. Moreover, be¬ 
cause of the large heat capacity of water and the ocean’s 
dynamics, the sea surface temperature (SST) anomalies 
vary in the scale of months. Calculating DI for a point in 
the central pacific (NIN03 area) for different values of r 
yields the results shown in Fig. [3j The point considered 
is the same as in Fig. |T] (a-c) and moreover, the panel [1] 
(c) is the same as panel [3] (a). 

For r = 1 Figj3] (a) shows the central Pacific influ¬ 
enced by (in blue) the eastern Pacific and influencing (in 
red), presumably through atmospheric teleconnections, 
the global tropics and the extratropical Pacific ocean. 
However, as r grows the number of significant connec- 


FIG. 4. (Color Online), (a), (b) Directionality of the signif¬ 
icant links (t = 1 month) of two nodes in the extratropics 
indicated with triangles: (a) in southern South America (de 
la Plata basin) and (b) in Labrador Sea. As the decorrelation 
time on the extratropics is very fast, no clear structures are 
seen. In order to show that it is not a problem of statistics, 
panels (c) and (d) display the significant MI values. 


tions decreases, suggesting that the time-scale of decor¬ 
relation of the SAT is less than 6 months. This is consis¬ 
tent with the persistence time scale of 3 to 6 months of 
observed sea surface temperature. 

The extratropical atmosphere shows larger internal 
variability than the tropics and the impact of the ex¬ 
tratropical SST on the atmosphere is much more limited 
than in the tropics. Thus, the variability of extratropi- 
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cal SAT is dominated by synoptic atmospheric dynamics 
and has time scales of a few days. Longer persistence 
time scales might appear in the extratropics if the re¬ 
gion is influenced by tropical SST. This motivates the 
use of a small value of r when considering extra-tropical 
variability. 

In Fig. 2] the DI and MI maps for two points in the 
extratropics are shown. Panels 0 ] (a,c) show links related 
to a point in southeastern South America, while panels 0 ] 
(b,d) show links related to a point in the Labrador sea, 
whose characteristics are linked to the North Atlantic 
Oscillation (NAO)p35|. The top panels show DI for r = 1 
month while the bottom panels show MI. 

Consistent with the previous description, the extra- 
tropical SAT show only some incoming links from the 
tropical region for r = 1 month. The point over the 
Labrador Sea seems to show also some outgoing links to 
the northeast, although there is no clear structure. As 
shown in Sec. IIII Cl a better identification of link direc¬ 
tionality is obtained by using daily data. 


2. Comparison between monthly and daily datasets 


(aj 



0 45 90 135 180 135 90 45 




0 45 90 135 180 135 90 45 


. 0.30 
^ 0.20 
— 0.10 
— 0.00 
—- 0.10 
1 - 0.20 
^- 0.30 



Ao.30 
® 0.20 
® 0.10 
— 0.00 
-0.10 


■-0.30 



FIG. 5. (Color Online). Comparison of monthly-averaged 
and daily-averaged datasets. The nodes considered are as in 
Fig. m Panels (a) and (c) display the DI measure (r = 1 
month) calculated from monthly data —same as Figs. [T| (c) 
and (f). Panels (b) and (d) display the DI measure (r = 30 
days) calculated from daily data. Results are consistent and 
the resolution using daily data is better. 


In order to obtain more temporal resolution daily data 
has been used. Figure [5] shows a comparison between 
DI for monthly data -panels [5] (a) and [5] (c)—and for 
daily data—panels [5] (b) and[5](d). Corresponding to the 
same points considered in Fig. [1] In order to adequately 
compare the datasets, r was chosen equal to one month 
in the monthly data and and 30 days in the daily data. 

The maps using monthly and daily data show similar 


features and no inconsistencies are found. Furthermore, 
the map constructed using daily data captures much bet¬ 
ter the local and remote dependencies and directional¬ 
ity of the links. Areas with significant links are better 
defined and some regions that are known to be influ¬ 
enced by equatorial Pacific SST, like the tropical north 
Atlantic [36[, clearly appear using daily data, but only 
very roughly using monthly data. Thus, the increase in 
temporal resolution improves the representation of the 
links related to tropical regions. 


C. Analysis of daily-averaged SAT anomalies 

1. Influence of the parameter r 
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FIG. 6. (Color Online). Influence of the parameter r when 
DI is computed from daily-averaged data. The DI value is 
plotted vs. r for the four links that connect a node in the 
Pacific ocean with (a) a node in Labrador sea, (b) a node in 
tropical north Atlantic (in the red area in Fig. [5] between the 
Caribbean and Africa), (c) a node in southern Pacific (the 
red spot in Fig. [5] near Antarctica) and (d) a node in the 
Indian ocean. In order to filter the noise the DI value was 
averaged over second neighbors in the grid. The shaded area 
(red online) indicates the DI value computed from 100 BS 
surrogates, as explained in section MM 


To analyze the influence of the parameter r in the DI 
measure when it is computed from daily data, four links 
connecting the node in the Pacific ocean are considered 
in Figs. 1, 5 with: 

(a) a node in Labrador sea, 

(b) a node in tropical north Atlantic (in the red area 
in Fig. [5] between the Caribbean and Africa), 

(c) a node in southern Pacific (the red spot in Fig. [5] 
near Antarctica) and 

(d) a node in the Indian ocean. 

Figure [6] displays, for each link, the DI value vs r. 

In panel (a), the directionality index that characterizes 
the influence of the equatorial Pacific to the Labrador 
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Sea shows a small increase for small r that persists up to 
about r = 10 days. This is the typical time scale associ¬ 
ated with the setup of the atmospheric anomalies forced 
by anomalies in tropical convection. Afterwards, the DI 
decreases exponentially-like, becoming non-significant at 
r = 30 days—as the DI value enters the shaded area. 
In the tropical north Atlantic—Fig. [6] (b)—the DI to 
the equatorial Pacific shows similar values as for the 
Labrador Sea case, for up to r = 10 days. For larger val¬ 
ues of r, however, the DI has significant values (approx¬ 
imately constant) up to about r = 60 days. Afterwards, 
it decreases becoming non-significant for r larger than 4 
months. The difference in behavior between the tropical 
north Atlantic and the Labrador Sea might be because 
the remote forcing from ENSO induces a clear regional 
response in the surface temperatures of the tropical At¬ 
lantic [36], which will add persistence to the remote sig¬ 
nal. On the other hand, in the Labrador Sea, the ocean 
does not respond and the large atmospheric variability 
obscures the signal from the equatorial Pacific. 


A similar behavior to the one seen in the tropical At¬ 
lantic is also found in the Southern ocean—Fig. [6] (c). 
There is a fast time scale for small values of r, but in 
this case the influence of the equatorial Pacific persists 
for r up to 80 days, with DI values significantly larger 
than in the Atlantic. The ENSO influence over the south 
Pacific is one of the most robust signals in the extratrop¬ 
ics, consequence of atmospheric teleconections associated 
with the Pacific South American pattern [37]. The time 
scale of about 3 months seen in the DI is likely associated 
with the time it takes to the surface ocean to respond to 
anomalous atmospheric fluxes and to the seasonal depen¬ 
dence of the atmospheric teleconnection pattern on the 
mean state of the extratropical atmosphere. 


The behavior of the influence of the equatorial Pacific 
onto the Indian ocean—Fig. [6] (d)—also shows a fast time 
scale of a few days, but in this case the largest value of the 
DI is seen for a r of about 60 days. Also, there are large 
values of DI for r values of more than 4 months. This 
suggests that the Indian ocean responds to the incoming 
ENSO signal in a time scale of about 2 months, through 
thermodynamic and dynamic coupling [38|, . 


Thus, the DI plots as a function of r for different ocean 
regions suggest that there is one fast time scale of about 
10 days associated to the setup of the atmospheric tele¬ 
conections associated with changes in tropical convec¬ 
tion. Moreover, in some regions the DI shows a second 
longer time scale of about 2 or 3 months associated with 
the response of the local ocean to the circulation anoma¬ 
lies forced from the tropics. Finally, longer time scales 
related to oceanic inertia also affect the DI value but 
longer datasets are needed in order to perform a robust 
estimation of their effects. 





FIG. 7. (Color Online). Effect of r in the southern extrat¬ 
ropics, when DI is computed from daily data. The same node 
in Fig. [4] (a) is considered, and the values of r are: (a) 1 day, 
(b) 3 days, (c) 7 days, and (d) 30 days. Small r capture wave 
trains propagating from west to east, while for r values longer 
than a week the decorrelation is higher and only the influence 
from the Pacific ocean persists. 



FIG. 8. (Color Online). As in Fig [ 7 | but for a node in the 
northern extratropics —the Labrador Sea, as in Fig. [4](b). 
The values of r are as in Fig [TJ The Labrador sea area is 
related to a source of atmospheric variability of the north 
Atlantic ocean that affects Europe. As in Fig. (7) smaller 
values of r can capture wave trains propagating from west 
to east, over Europe, while values of r longer than a week 
loose this effect and only the influence from the Pacific ocean 
persists. 


2. Influence of r in the extra-tropics 

The improvement in characterizing directionality us¬ 
ing daily data is even larger in the extratropics. As men¬ 
tioned above, the extratropical SAT is strongly depen¬ 
dent on synoptic scale perturbations (a time scale of a 
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few days). Thus, the use of daily data should allow to 
uncover these relationships and investigate the direction 
of the links as the lag increases. In order to do so, the 
point in southeastern South America shown in Fig. 0 ] 
is considered and the directionality network for several 
values of r ranging from a day to one month (Fig. [7j) 
is constructed. For synoptic time scales of a few days 
the methodology uncovers the existence of a wave train 
connected to southeastern South America propagating 
with a southwest-northeast direction. Moreover, there 
is a clear separation line between regions with incoming 
and outgoing links. This configuration is characteristic 
of the progression of a front through the reference point 
and does not imply that the SAT over the reference point 
influences the region to the northeast but it only happens 
to be in the path of the perturbation. As the lag time in¬ 
creases, the extratropical wave train associated with syn¬ 
optic time scales fades and only the points in the tropics 
remain, consistent with an influence of the equatorial Pa¬ 
cific on the region on longer time scales, perhaps related 
to ENSO. 

A similar behavior is seen taking as reference point 
the SAT over the Labrador Sea (Fig. [8]). For small val¬ 
ues of r the progression of a front is clearly detected 
using this procedure: given the mean westerly winds at 
these latitudes, the front moves from west to east and 
is clearly marked as the boundary between the incoming 
and outgoing links. It is also seen for r = 3 suggesting 
that in about three days the front reaches the Mediter¬ 
ranean region affecting temperatures there. Again, as r 
increases mainly the tropical links remain. However, even 
for r = 30 there is a well defined region of outgoing links 
that remain over the Labrador Sea, suggesting that the 
SAT in the region may have relatively long time scales 
of variability, perhaps related to the North Atlantic Os¬ 
cillation. Another possibility is the influence of the local 
ocean that increases the persistence of atmospheric tem¬ 
perature anomalies through thermodynamic coupling, as 
shown by Barsugli et. al. [40| . 


3. Influence of r in the tropics 

Figure [9] presents the effect of r on daily data in the 
tropics. Panel [9] (a) shows the point considered in the 
Pacific ocean for r = 1 day. Contrarily to e.g. figure [7J 
where front propagation shows a big response in the trop¬ 
ics due to the existence of wave trains in the extratropics, 
in this case the connections are weak. The opposite case, 
for r > 30 days—Panels [9] (b,c)—are consistent with in 
figure [6] as the structures are robust for a wide range of 
r. Notice in Panel [9] (c) the presence of a second blue 
area near the coast of Chile, This is in agreement with 
the easterly trades advection over the Pacific ocean. In 
order to complement the information provided in panels 
[9] (a,b,c) and also in Fig. [6j in the supplementary in¬ 
formation a movie was included 0 (Multimedia view), 
which presents the links to the node in the Pacific ocean, 


for r = 1 day to r = 180 days. Panels [9] (d,e,f) show 
maps for a point in the Indian ocean, for values of r = 1, 
45, and 90 days respectively where similar results are ob¬ 
tained. 


Directionality on the tropical Pacific Ocean 

In Fig. [10] the DI for r = 30 days has been plotted for 
points covering all the equatorial Pacific. They begin at 
95.0° W -panel ITOl (a) - and end at 125.0° E -panel HOl 
(d). Clearly, the influence the Pacific ocean exerts is al¬ 
most global, over tropics and extratropics, in agreement 
with previous studies (e.g. in [Hf). The DI allows to 
show that —even if there are feedbacks and the Pacific is 
affected by extratropical perturbations and other ocean 
basins (e.g. the tropical Atlantic)— the influence is effec¬ 
tively from the Pacific to the rest of the world. Moreover, 
the maps show that the largest influence is done by the 
equatorial Pacific close to the dateline. This is clear in 
the extratropical atmosphere, as well as in the tropical 
north Atlantic. On the other hand, the connection to 
the Indian ocean and south Atlantic is not so sensitive 
of the point considered over the equatorial Pacific. As 
the reference point moves further west from the date¬ 
line the influence decreases substantially, only remaining 
weak connections to the tropical north Atlantic and In¬ 
dian oceans. The methodology can thus be applied to 
find the best region to construct an index that describes 
the Pacific influence over the area where climate anoma¬ 
lies are studied. 

Notice that all maps show a blue tongue of incoming 
links to the east of the point considered; it it is seen first 
in panel 00] (a) and extends westward until covering the 
whole Pacific ocean in 00 ] (d). This feature is related to 
the existence of the equatorial cold tongue and the fact 
that easterly trades blow over the equator thus advecting 
air from the east to the west of the point. 


IV. CONCLUSIONS 

The climatological relevance of directed climate net¬ 
works constructed via information-theoretic tools has 
been shown. The presence of significant links was in¬ 
ferred using mutual information while the direction of 
those links was inferred using the directionality index. 
Using monthly-averaged SAT data the results from our 
previous work |33j have been recovered and the analysis 
extended to infer the net direction of propagation of the 
information. The inference method was tested against 
the value of the parameter r, that represents the time 
required for the information to travel from one node to 
another. It was found that by adequately tuning r the 
network connectivity varied revealing the various time- 
scales of atmospheric processes; too short values of r 
failed to capture several long-range links, while for too 
large values of r (above the length of correlations in the 
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FIG. 9. (Multimedia view) (Color Online). Effect of r using daily data for a node in the tropics. In panels (a-c) the node is 
in the Pacific ocean, and (d-f), in the Indian ocean. The values of r are: (a,d) 1 day, (b,e) 45 days, and (c,f) 90 days. The 
maps for r = 30 were shown in Fig. [3] In contrast to the extra-tropics, where the propagation of waves in the scale of days is 
dominant, and are clearly seen in the maps for time scales of a few days, in the tropics the variability is longer as the ocean 
adds persistence to the nearly barotropic atmosphere. For higher values of r only long lasting phenomena are observed, related 
to the strongest variability pattern in this area, which is ENSO. In the supplementary information a movie showing the links 
to the node of the fist row, for r = 1 to r = 180 is presented [ill ]. 


data) the connectivity of the network decreased drasti¬ 
cally. 

In addition, when considering daily-averaged SAT 
data, the analysis revealed variability patterns consis¬ 
tent with known features of the global climate dynam¬ 
ics. In the extratropics the long time average synoptic 
weather was correctly inferred: as specific examples two 
geographical regions in different hemispheres—one node 
in de la Plata basin and another in the Labrador Sea— 
were considered, and the link direction revealing wave 
trains propagating from west to east, in both hemispheres 
was shown. 

As a future work, it would be interesting to apply this 
methodology to analyze how the climate network changes 


in the different seasons. Another interesting issue is the 
role of permutation entropy. Computing the direction¬ 
ality index using the ordinal patterns approach nisi, 
which would allow constructing networks revealing atmo¬ 
spheric processes with short time-scales (of days-months) 
or with long time-scales (few years). These studies are in 
progress and will be reported elsewhere. 
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